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o 
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■ Abstract. We study an integrable case of n-particle Toda lattice: open chain with boundary 

terms containing 4 parameters. For this model we construct a Backlund transformation 
and prove its basic properties: canonicity. commutativity and spcctrality. The Backlund 
transformation can be also viewed as a discretized time dynamics. Two Lax matrices are 
used: of order 2 and of order 2n + 2, which are mutually dual, sharing the same spectral 
curve. 
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In the present paper we study the Hamiltonian system of n one-dimensional particles with 



1 Introduction 

(N 
> 

ly-j \ coordinates Xj and canonical momenta Xj, j = 1, . . . , re: 

o\ ' 

{Xj,X k } = {xj,x k } = 0, {Xj,x k } = 5j k , (1.1) 
' characterized by the Hamiltonian 

o 



i 1 1 

H = ^T -X* + J2 e Xj+1 ~ Xj + «ie Xl + -/3ie 2xi + a n e~ Xn + -(3 n e~ 2x " (1.2) 



n-l 

>< 

containing 4 arbitrary parameters: a±, Pi, a n , f3 n . 

The model was missing from the early lists of integrable cases of the Toda lattice [H [2] based 
on Dynkin diagrams for simple affine Lie algebras. Its integrability was proved first in [3j[U[5]. 
As for the more recent classifications, in [6] the model is enlisted as the case (i). In [TJ [8] 
particular cases of the Hamiltonian (jl.2p are assigned to the Cn^ case with 'Morse terms'. For 



brevity, we refer to the model as 'BC-Toda lattice' emphasising the fact that each boundary 
term is a linear combination of the term ~ a corresponding to the root system B and of the 
term ~ f3 corresponding to the root system C, see [U EJ [TJ E] . 

In section [2] we review briefly the known facts about the integrability of the model using 
the approach developed in [31 H] and based on the Lax matrix L{u) of order 2 and the corre- 
sponding quadratic r-matrix algebra. In particular, we construct explicitly a generating function 
of the complete set of commuting Hamiltonians Hj (j = 1, . . . ,n) which includes the physical 
Hamiltonian H (fl~2|) . 



*This paper is a contribution to the Vadim Kuznetsov Memorial Issue 'Integrable Systems and Related Topics'. 
The full collection is available at |http://www.emis.de/journals/SIGMA/kuznetsov.html| 
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In Section [3] we describe the main result of our paper: construction of a Backhand transfor- 
mation (BT) for our model as a one-parametric family of maps B\ : (Xx) \— > (Yy) from the 
variables (Xx) to the variables (Yy). We construct the BT choosing an appropriate gauge (or 
Darboux) transformation of the local Lax matrices. In Section [U adopting the Hamiltonian 
point of view developed in [SCEO], we prove the basic properties of the BT: 

1. Preservation of the commuting Hamiltonians B\ : Hj(X,x) t— > Hj(Y,y). 

2. Canonicity: preservation of the Poisson bracket (jl.ip . 

3. Commutativity: B\ 1 o B\ 2 = B\ 2 o B\ 1 . 

4. Spectrality: the fact that the graph of the BT is a Lagrangian manifold on which the 
2-form 

n 

n = J2{dXj A dxj - dYj A dyj) - din A A dA (1.3) 

3=1 

vanishes. Here A is an eigenvalue of the matrix L(\). In other words, the parameter A of 
the BT and its exponentiated canonical conjugate A lie on the spectral curve of L(u): 

det(A-L(A)) = 0. (1.4) 

We also prove the following expansion of B\ in A -1 

B x :f^f-2\~ 1 {HJ} + 0(\~ 2 ), A^oo. (1.5) 

which allows to interpret the BT as a discrete time dynamics approximating the continuous-time 
dynamics generated by the Hamiltonian (|1.2|) . 

In Section [5] we construct for our system an alternative Lax matrix C(v). The new Lax 
matrix of order 2n + 2 is dual to the matrix L(u) of order 2 in the sense that they share the 
same spectral curve with the parameters u and v having been swapped: 

det(v -L(u)) = (-l) n+1 vdet(u- C(v)). (1.6) 

In the same section we provide an interpretation of the BT in terms of the 'big' Lax mat- 
rix C(v) and establish a remarkable factorization formula for A 2 — C 2 (v). 

The concluding Section [6] contains a summary and a discussion. All the technical proofs and 
tedious calculations are removed to the Appendices. 

2 Integrability of the model 

In demonstrating the integrability of the model we follow the approach to the integrable chains 
with boundary conditions developed in [31 H] and use the notation of [9l [10] . 
The Lax matrix L(u) for the BC-Toda lattice is constructed as the product 

L(u) = J<'_(u)T t (-'u)J<' + (u)T(u) (2.1) 

of the following matrices (T* stands for the matrix transposition). 
The monodromy matrix T(u) is itself the product 

T(u) = £ n ( u ) ■ ■ ■ h(u) (2.2) 

of the local Lax matrices 

l 3 (u) = £(u; X 3 , Xj ) = ( u + , (2.3) 
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each containing only the variables Xj, xj describing a single particle. Note that trT(tt) is the 
generating function for the Hamiltonians of the periodic Toda lattice. 

The matrices K±{u) containing the information about the boundary interactions are defined 

as [am 

K-(u) = ( U ~ a A t K + (u) = ( U -*»). (2.4) 

The significance of the Lax matrix L(u) is that its spectrum is invariant under the dynamics 
generated by the Hamiltonian (|1-2|) . the corresponding equations of motion dG/dt = G = {H, G} 
for an observable G being 

± j = x h J = l,...,n (2.5) 

and 

Xj = tfi+i-*' - e 3 * -3 *- 1 , j = 2, . . . ,n - 1, (2.6a) 

Xi = e* 2 "* 1 - ai e Xl - foe 2 * 1 , (2.6b) 

X n = -e*™-*™- 1 + a n e- x " + /? n e- 2x " . (2.6c) 

To prove the invariance of the spectrum of L(u) we introduce the matrices Aj(u) 

M^) = [_~~l- 1 6 q)' i = 2,...,n-l, (2.7) 

a ( \ ( ~ u qXi \ a ( \ ( ~ u ot n + (3 n e- Xn \ , > 

Mu) =[-ai-Pie*i oj' A ^^ = {-e-* n J' (2 " 8) 

which satisfy the easily verified identities 

lj = A j+1 £j - ijAj , j = 1, . . . , n, (2.9) 
-£+ = = K+^ n+ i(n) + A t n+1 (-u)K+, (2.10a) 
/<■_ = = Ax{u)K- + if_^(-u). (2.10b) 



From (12. 2j) and (12, 9j) it follows immediately that 

= A n+1 (u)T(u) - T{u)Ai{u). (2.11) 

Then, using (12. 1|) and (I2.10|) . we obtain the equality 

L(«) = [Ai(u),L(u)] (2.12) 

implying that the spectrum of L{u) is preserved by the dynamics. 

There are only two spectral invariants of a 2 x 2 matrix: the trace and the determinant. 
From (|2.3[) it follows that det i{u) = 1 and, respectively, det T(u) = 1, so, by (12. If) . the determi- 
nant of L(u) 

d{u) =detL(n) = det K- (u) det K+ (it) = (a? + Piu 2 )(a 2 n + /? n u 2 ) (2.13) 
contains no dynamical variables Xx. The trace 

t(u) = ixL(u) = trK_(n)T*(-u)K + (n)T(n), (2.14) 
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however, does contain dynamical variables and therefore can be used as a generating function of 
the integrals of motion, which can be chosen as the coefficients of the polynomial t{u) of degree 
2n + 2 in u. Note that t(—u) = t(u) due to the symmetry 

K^i-u) = -K±{u). (2.15) 

The leading coefficient of t(u) at u 2n+2 is a constant (— l) n . Same is true for its free term 

i(0) = tr K + (0)K-(0) = -2a n ai (2.16) 

due to the identity 

MK ± (0)M t = detM ■ K±(0), (2.17) 

which holds for any matrix M. 

We are left then with n nontrivial coefficients Hj 

n 

t( u ) = {-l) n u 2n+2 - 2a n ai + ^ H jU 2j (2.18) 

j'=i 

which are integrals of motion Hj = since t(u) = due to (|2.12p . 

The conserved quantities Hj are obviously polynomial in X, e ±x . Their independence can 

easily be established by setting e ±x = in (|2.3[) and analysing the resulting polynomials in X. 
It is also easy to verify that the physical Hamiltonian (|1.2p is expressed as 

(_l)n+l 

H = ^ H n . (2.19) 

The quantities Hj are also in involution 

{Hj,H k } = (2.20) 

with respect to the Poisson bracket (11. ip . Together with the independence of Hj, it constitutes 
the Liouville integrability of our system. 

The commutativity (|2.20j) of Hj or, equivalently, of t(u) 

{t(u 1 ),t(u 2 )} = (2.21) 

is proved in the standard way using the r-matrix technique [31 S] ■ 
Let 1 be the unit matrix of order 2 and for any matrix L define 

l 2 

L = L<g)l, L = 1<S>L. (2.22) 
We have then the quadratic Poisson brackets [TOj, [11] 

{}{ui) ,£{u 2 )} = [r(«i-«2),J(«i)?(«2)], (2.23) 
and, as a consequence, 

{T{ui), T(u 2 )} = [r(«i - u 2 ), T( Ul )T(u 2 )l (2.24) 
with the r-matrix 

r(u) = -, (2.25) 
u 

where V is the permutation matrix Va ®b = b® a. 
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Let 

r(u) = r h (u) = r t2 (u), (2.26) 

t\ and t 2 being, respectively, transposition with respect to the first and second component of 
the tensor product C 2 ® C 2 . 

Then for both T(u) = T(u)K_(u)T t (—u) and T(u) = T f (— u)K + (u)T(u) we obtain the same 
Poisson algebra [31 H] 

12 12 12 

{T(ui),T(u 2 )} = r(u x - u 2 )T(u 1 )T(u 2 ) - T(u x )T(u 2 )r(ui - u 2 ) 

1 2 2 1 

- T( Ul )r( Ul + u 2 )T(u 2 ) + T(u 2 )r( Ul + u 2 )T{ Ul ), (2.27) 
which ensures the commutativity (|2.21|) of t(u). 

3 Describing Backlund transformation 

In this section we shall construct a Backlund transformation (BT) for our model. We shall stay 
in the framework of the Hamiltonian approach proposed in [9] and follow closely our previous 
treatment of the periodic Toda lattice [9j [10], with the necessary modifications taking into 
account the boundary conditions. 

We are looking thus for a one-parametric family of maps B\ : (Xx) i— > (Yy) from the 
variables (Xx) to the variables (Yy) characterised by the properties enlisted in the Introduction: 
Invariance of Hamiltonians, Canonicity, Commutativity and Spectrality. 

The invariance of the commuting Hamiltonians Hj, or of their generating polynomial t(u) = 
trL(u) will be ensured if we find an invertible matrix M\(u, A) intertwining the matrices L(u) 
depending on the variables Xx and Yy: 

Mi(u,A)L(u; Y,y) = L(u;X,x)M 1 (u, X). (3.1) 
To find M\ (u, A) let us look for a gauge transformation 

M j+1 (u,X)£(u;Y j ,y j )=£(u;X j ,x j )M j (u,X), j = l,...,n, (3.2) 
implying that detMj does not depend on j. From (|3.2[) and (|2.2p we obtain 

M n+1 (u,X)T(u;Y,y) = T(u;X,x)M 1 (u,X). (3.3) 
Let J be the the standard skew-symmetric matrix of order 2 

J=( 1 n V J' = ~J, J 2 = ~l, (3.4) 



x 1 , 

and define the antipode M a as 

M a = -JMJ (3.5) 
for any matrix M of order 2. It is easy to see that 

M*M° = M a M l = det M. (3.6) 



Transposing (13. 3ft and using f)3.6f) together with the the fact that det Mj is independent of j 
we obtain the relation 

T t (-u;X,x)M- +1 (-u,X) =M?(-u,X)T t (-u;Y,y). (3.7) 
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We shall be able to obtain (|3.ip if we impose two additional relations 

K-(u)Mf(-u,\) = M 1 (u,X)K^(u), (3.8a) 
K + (u)M n+1 (u, A) = M% +1 (-u, X)K+(u). (3.8b) 



Then, starting with the right-hand side L(u; X, x)M±(u, A) of (|3.ip and using (|2.ip and ()3.3p 
we obtain 

L(u;X,x)Mi(u,A) = K-{u)T t (-u;X,x)K + (u)T(u;X,x)M 1 (u,X) 

= K-(u)T\- U ;X,x)K + ( U )M n+1 (u,X)T(u;Y,y) (3.9) 

Using then (ETBbl) to move M n+ i(ti, A) through K + (u), then using (|3.7p and finally (|3.8ap we 
get, step by step, 

L(u; X, x)M 1 (u, A) = K- (u)T* (-«; X, x)M£ +1 (-n, A) if+ (u)T(u; Y, y) 
= ET„( U )Af 1 t (- U , A)r*(-u; F y)K + (u)T(u; F, y) 
= Mi (u, X)K_ {u)T\-u- F, y)K + (u)T(u; Y, y) 

= M 1 {u,X)L(u;Y,y) (3.10) 



arriving finally at (13. ip . 

We have thus to find a set of matrices Mj(u, A), j = l,...,n + 1 compatible with the 
conditions (13. 2p and (|3.8p . A quick calculation shows that the so called DST-ansatz for Mj used 
in [9l [10] for the periodic Toda lattice contradicts the conditions (|3.8p . 

The philosophy advocated in [10] requires that the ansatz for the gauge matrix Mj{u) be 
chosen in the form of a Lax matrix satisfying the r-matrix Poisson bracket (|2.23p with the same 
r-matrix (|2.25p as the Lax operator £(u). It was shown in [10] that the so-called DST-ansatz 

Mf T («,A)=( n " A + S ^ (3.11) 

serves well for the periodic Toda case. The above ansatz is however not compatible with the 
boundary conditions (|3.8p and we have to use a more complicated ansatz for Mj in the form of 
the Lax matrix for the isotropic Heisenberg magnet (XXX-model): 

Mj(u, A) = ( U " A + SjSj s2 J Sj 7 2A % ) , det Mj (u, A) = A 2 - u 2 . (3.12) 
V —U — X + SjbjJ 

The same gauge transformation was used in [12] for constructing a Q-operator for the quan- 
tum XXX-magnet. 

Substituting (13. 12H into (|3.2p we obtain the relations 



Xj = -X + sjV' + Sj+ie - ^ , (3.13a) 

Yj = X - srV^ - Sj +1 e~ y i , (3.13b) 

Sj = 2XSJ 1 - sj 2 e x i - sj 2 e y i , (3.13c) 

S j . | e ' ■ • c . (3.13d) 

for j = 1, . . . , n, and from (|3.8p . respectively, 

2(ai + ftAsi) 2(As ra+ i - a n ) 
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Eliminating the variables Sj, we arrive to the equations defining the BT (j = 1, . . . ,n): 

Xj = -A + sj V J + s j+ ie~ Xj , (3.15a) 
Yj = A - sjV* - s j+1 e~ y i . (3.15b) 



The variables Sj, j = 1, . . . ,n + 1 in (|3. 15 j) are implicitly defined as functions of x, y and A 
from the quadratic equations 

(e"^- 1 + e-W- 1 )^ - 2Asj + (e^ + e») = 0, j = 2, . . . , n (3.16a) 
(2«i + fte 11 + P x e yi )s\ - 2\ Sl + (e xi + e^ 1 ) = 0, (3.16b) 
(<sT Xn + e- yn )s 2 n+1 - 2As n+ i + (2a n + (3 n e~ Xn + p n e~ yn ) = 0. (3.16c) 



Like in the periodic case [91 [TO], the BT map B\ : (Xx) i— > (Yy) is described implicitly by 
the equations (|3.15p . Unlike the periodic case, we have extra variables Sj. It is more convenient 
not to express Sj from equations (|3.16p and to substitute them into f)3. 15|) but rather define the 
BT by the whole set of equations f)3. 15H and (I3.16|) . 

Equations f)3. 15|) and p,16p are algebraic equations and therefore define (Yy) as multivalued 
functions of (Xx), which is a common situation with integrable maps |13j . 

In this paper, to avoid the complications of the real algebraic geometry we allow all our 
variables to be complex. 



4 Properties of the Backlund transformation 

Having defined the map B\ : (Xx) t— > (Yy) in the previous section, we proceed to establish its 
properties from the list given in the Introduction. 

4.1 Preservation of Hamiltonians 

The equality Hj(X, x) = Hj(Y, y) V A, or, equivalently, t(u; X, x) = t(u; Y, y) holds by construc- 
tion, being a direct consequence of (13. ip . 

4.2 Canonicity 

The canonicity of the BT means that the variables Y(X, x; A) and y(X, x; A) have the same 
canonical Poisson brackets (jl.ip as (Xx). An equivalent formulation can be given in terms of 
symplectic spaces and Lagrangian manifolds. Consider the 4n-dimensional symplectic space V/^ n 
with coordinates XxYy and symplectic 2-form 

n 
3=1 

Equations ()3. 15j) and (|3.16p define a 2n-dimensional submanifold T2 n C V^ n which can be 
considered as the graph Y = Y(X, x; A), y = y(X, x; A) of the BT (the parameter A is assumed 
here to be a constant). The canonicity of the BT is then equivalent to the fact that the manifold 
T2n is Lagrangian, meaning that: (a) it is isotropic, that is nullifies the form £l^ n 

ft 4 nlr 2 „ = 0, (4.2) 



and (b) it has maximal possible dimension for an isotropic manifold: diml^n = idim!4n- 
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One way of proving the canonicity is to present explicitly the generating function &\(y; x) of 
the canonical transformation, such that 

The required function is given by the expression 

n 

$\{y,x) = ^2f\(yj,sj+i;xj,sj) + ^ 0) («i) + ^ n+1) (s„+i), (4.4) 

3=1 

where 

f\(yj,Sj+r, Xj,sj) = X(2lnsj - xj - yj ) + sj l (e x i + - s j+ i(e - ^ + e"%), (4.5a) 

2ai 



</?f(si) = -Aln(l + /?is?) - ^arctan( v / Asi), (4.5b) 

VPl 

^l n+1) ( s ^+l) = Aln (Ai + 4+i) - -^arctan ( ^7= ) , (4.5c) 

V Pn \ V Pn J 

and Sj(x,y; A) are defined implicitly through (|3. 16j) . 

Equalities (|4.3jl can be verified by a direct, though tedious, computation. Another, more 
elegant, way is to use the argument from [lOj based on imposing a set of constraints in an 
extended phase space, see Appendix O 

4.3 Commutativity 

The commutativity B\ 1 oB\ 2 = B\ 2 oB\ L of the BT follows from the preservation of the complete 
set of Hamiltonians and the canonicity by the standard argument [21 ID] based on Veselov's 
theorem [13] about the action-angle representation of integrable maps. 



4.4 Spectrality 

The spectrality property formulated first in [9] generalises the canonicity by allowing the pa- 
rameter A of the BT to be a dynamical variable like x and y. 

Let us extend the symplectic space V^ n from section 14.21 to a (4n + 2)-dimensional space 
^4n+2 by adding two more coordinates A, \i and defining the extension ^l^ n+ 2 of symplectic 
form Q,4 n ([4.1 p as 

n 

= ^An - d/i A dA = (dXj A dxj - dYj A dyj) — d/i A dA. (4.6) 
i=i 

Define the extended graph I^n+i of the BT by equations ([3.15P and a new equation 
d 

^ = -^x(y,x). (4.7) 

The 2- form 04 n+ 2 obviously vanishes on I^n+i, and the manifold T2 n +i is lagrangian. 
An amazing fact is that e^ 1 is proportional to an eigenvalue of the matrix L(X), see f)1.4|) . In 
fact, the two eigenvalues of L(\) can be found explicitly to be 

A=(al + (3 n \*)p^fl(-sfe^), (4.8a) 

Pn + s n+l j = i 



Backlund Transformation for the BC-Type Toda Lattice 



9 



A = (a? + AA 2 )^-^ IIH 6 ^)' ( 48b ) 
i -+- Pi^i . =1 

see Appendix |B] for the proof. 

Having the explicit formulae (|4,8ap for A and (|4,4p for <I>^(y;rc) one can easily verify that 

A = (-ir(o 2 ; + /3 n A 2 )e^. (4.9) 
4.5 Backlund transformation as discrete time dynamics 

One of applications of a BT is that it might provide a discrete-time approximation of a continuous- 
time integrable system [14] 115] . Indeed, iterations of the canonical map B\ generate a discrete 
time dynamics. Furthermore, if we find a point A = Ao that (a) the map B\ becomes the iden- 
tity map, and (b) in a neighbourhood of Ao the infinitesimal map B\ Q+e ~ e{H, •} reproduces 
the Hamiltonian flow with the Hamiltonian (|1.2j) . we can claim that B\ is a discrete time ap- 
proximation of the BC-Toda lattice. An attractive feature of this approximation is that, unlike 
some others [14] , the discrete-time system and the continuous-time one share the same integrals 
of motion. 

In our case Ao = oo. Letting e = A -1 and assuming the ansatz 

Vi = Xj + 0{e), j = l,...,n (4.10) 
we obtain from (|3. 16a.]) and (|3.16bp the expansion 

Sj = ee Xj + 0(e 2 ), j = 1, . . . , n (4.11a) 

and from (|3.16cj) the expansion 

s n+ i = e{a n + (3 n e- x -) + 0(e 2 ). (4.11b) 
Substituting then expansions (pTTU]) into equation (|3.13dD we obtain 

Sj = 2e~ x ^ 1 +0(e), j = 2, ...,n+l (4.12a) 

and substituting expansion (|4.1ip for s± into formula (|3.14j) for S\ we obtain 

S 1 = (2a 1 + (3 1 e x ')+0(e). (4.12b) 

Then from (|3.12p we have 

-eMj = l + e(ul + 2Aj) +0(e 2 ), j = l,...,n + l, (4.13) 

where Aj coincides with the matrix (given by (|2.7p and (|2.8|) ) which describes the continuous- 
time dynamics of the Lax matrix. From (|3.2p we obtain then 



£(u; Yj,yj) = £(u; Xj, Xj ) - 2e(A j+1 e(u; Xj, Xj ) - £(u; X^x^Aj) + 0(e 2 ), (4.14) 
for j = 1, . . . , n + 1. Comparing the result to (|2.9p we get the expansion (|1.5p . 
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5 Dual Lax matrix 



Many integrable systems possess a pair of Lax matrices sharing the same spectral curve with 
the parameters u and v swapped like in (jl.6p . see [16] for a list of examples and a discussion. 
In particular, the periodic n-particle Toda lattice has two Lax matrices: the 'small' one, of 
order 2 [TT], and the 'big' one, of order n |17j . For various degenerate cases of the BC-Toda 
lattice 'big' Lax matrices are also known [2] \7\ [8| [TT] . 

In this section we present a new Lax matrix of order 2n + 2 for the most general, 4-parametric 
BC-Toda lattice. Here we describe the result, removing the detailed derivation to Appendix O 

Let Ejk be the square matrix of order 2n + 2 with the only nonzero entry (Ej^jk = 1. The 
Lax matrix C(v) is then described for the generic case n > 3 as 



C(v) = E CjkEjk 

j,k=l 

n n 

= E ' ■ + E(-*y% + E u+i) 

=2 j=l 
n—1 n 

^g^j+i x i E2n+2-j,2n+l-j + E ip^j ^2n+2- j,2n+2- j ~ E2 n +2-j,2n+3-j) 
3=1 3=1 

a n e~ Xn + ye -2 *"^ (£ n +l,n - E n+2 ,n+l) 

n -x„-i ^E n+3 n _ E n+ 2,n-l) — E n+ i yU+ 2 
(a x e^ + Y^ 1 ) ( fi 2Wi+""^l,2n+2) 



+ ^e X ^(E 2 ,2n + l - E^n) ~ vE 2n+ 2,l 
2v 



(5.1) 



and consists of a bulk 'Jacobian' strip (the main diagonal and two adjacent diagonals) which 
reproduces the Lax matrix for the open Toda lattice together with boundary blocks containing 
parameters ot\f3\OL n (3 n . We do not consider here the special case of small dimensions n = 1, 2 
when the two boundary blocks interfere with each other and the structure of the Lax matrices 
becomes more complicated 

To help visualise the matrix C(v) we present an illustration for the case n = 3, using the 
shorthand notation = e Xj , rjj = e yj : 








-1 



C(v) 



(- 


-x 1 


1 




«2 
fi 


-x 2 







«3 



















f>3 
2«2«3 


















V 


— V 















01 1 t 





1 











0±t t 


-x 3 


1 
















-1 













"3 03 
€3 2£f 


X 3 


-1 





2€ 2 C3 





£3 

52 


x 2 


-1 











(2 

Si 



















(5.2) 



The matrix C(v) possesses the symmetry 
£(v) = -C v C\v)C-\ 



(5.3) 
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where 



2n+l 

C v = —vE2n+2,2n+2 + E Ej,2n+2-j 

3=1 












1 








. 


. 1 











1 . 


. 








1 


. 


. 











. 


. 





-v 



(5.4) 



(note that C~ l = C v -i). 

The matrix C(v) shares the same spectral curve with the 'small' Lax operator L{u) satisfying 
the determinantal identity (|1.6jl and thus generates the same commuting Hamiltonians Hj. 

The Lax matrix C(v) of order 2n + 2 seems to be new. When one or more of the constants 
a\(3ia n (3 n vanish it degenerates (with a drop of dimension) into known Lax matrices for simple 
affine Lie algebras [21 E El U] - For the general 4-parametric case a Lie-algebraic interpretation 
of C(v) is still unknown. In particular, it is an interesting question whether C(v) satisfies a kind 
of r-matrices Poisson algebra. 

Inozemtsev [5] presented a different Lax matrix for the BC-Toda lattice, of order 2n instead 
of 2n + 2 and with a more complicated dependence on the spectral parameter. The relation of 
these two Lax matrices is yet to be investigated. 

For the dynamics (|2.5|) , f|2.6j) we have an analog of the Lax equation (|2.12p : 



C{v) = [A(v),C(v)] 
with A(v) defined as 



(5.5) 



3=1 



J.) 



E 



3,3+1 



XjE 2n+2 _j 



ft 



,2n+2-j + E2n+2-j,2n+3-j) + ^n+l.n+l 
Pn 



^e 2 * 1 {E 2n+2M+1 + v- l E 1)2n+2 ) + t^ e ~^n (E n+1 , n - E n+2 , n ) (5.6) 



+ vE 2n+2 ^ 

and satisfying 

A(v)C v + C v A\v) = 0. 

The analog of the formula (|3.1|) for the Backlund transformation is 

M(v,X)C(v;Y,y) = £{v; X, x)M(v, A), 
M{v, X)C(v; X, x) = £{v; Y, y)M{v, A), 

where M{v) is given by 



j,k=l 

n-1 

+ £ 

3=1 

a 



3=2 



Vj-i 



3=1 



Sj+1 Cj 



n, 



E 33 + E 3,j+1 



(5.7) 

(5.8a) 
(5.8b) 

(5.9) 



Vj+i 
6 



3 



2n+2-j,2n+l 



+ 



3=1 

[E n +l,n — E n+2 ^ n+ ij + 



S 3+l _ Hi 



0n 



(Pis \ 
ttlflH Y~ j {E 2n +2,2n+l+ V ~ 1 E ian+2) + 



2£n£ri— 1 

2w 



E2n+2-j,2n+2-j ~ -£'2n+2-j,2n+3-j 



2,2n+l" 



2n+2,l> 
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(using again the notation £j = e Xj , rjj = e yj ) and M. (v) is defined as 

Mty^CvM^C- 1 . (5.10) 

One of common ways to obtain a Backlund transformation is from factorising a Lax matrix 
in two different ways, see [18] for Toda lattices and [13] for other integrable models. For our 
model we also have a remarkable factorisation, only instead of Civ) we have to take its square: 

A 2 - C 2 {v; X, x) = M(v, X)M(v, A), (5.11a) 
A 2 -C 2 (v;Y,y) = M(v, X)M(v, A). (5.11b) 



6 Discussion 

The method for constructing a Backlund transformation presented in this paper seems to be 
quite general and applicable as well to other integrable sZ(2)-type chains with the boundary 
conditions treatable within the framework developed in [31 H] . 

There is little doubt that a similar BT can be constructed for the D-type Toda lattice and 
a more general Inozemtsev's Toda lattice [5] with the boundary terms like 

ai h a n b n 

sinh 2 sinh 2 x\ sinh 2 ^ sinh 2 x n 

since those, as shown in [20], can also be described in the formalism based on the boundary K 
matrices (|2.ip and the Poisson algebra (|2.27j) . 

The 'big' Lax matrix C(v) still awaits a proper Lie-algebraic interpretation. Obtaining a BT 
from the factorisation of A 2 — £ 2 like in (|5,lip might prove to be useful for other integrable 
systems related to classical Lie algebras. 

It is well known that the quantum analog of a BT is the so-called Q-operator [21], see also [9]. 
Examples of Q-operators for quantum integrable chains with a boundary have been constructed 
recently for the XXX magnet [12] and for the Toda lattices of B, C and D types [22]. Our 
results for the BC-Toda lattice agree with those of [22], the generating function of the BT being 
a classical limit of the kernel of the Q-operator. Hopefully, our results will help to construct the 
Q-operator for the general 4-parametric quantum BC-Toda lattice. 



A Proof of canonicity 

Here we adapt to the BC-Toda case the argument from |10] developed originally for the periodic 
case. The trick is to obtain the graph T2 n of the BT as a projection of another manifold in 
a bigger symplectic space, the mentioned manifold being Lagrangian for trivial reason. 

Consider the 8-dimensional symplectic space Ws with coordinates XxYySsTt and the sym- 
plectic form 

uj 8 = dX A dx + dS A ds - dY A dy - dT A dt. (A.l) 
The matrix relation 



M(u, A; T, t)£(u; Y, y) = l(u\ X, x)M(u, A; S, s) (A.2) 
is equivalent to 4 relations 

X = -X + s~ l e x + te~ x , (A.3a) 
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Y = A - s _ V - te~ v , (A.3b) 
S = 2A.S" 1 - s~V - s~ 2 e y , (A.3c) 
T = e~ x + e"", (A.3d) 

defining a 4-dimensional submanifold ^4 C W%. The fact that Q4 is Lagrangian, that is ^s\g 4 = 0, 
is proved by presenting explicitly the generating function 

f x (y, t; x, s ) = X(2lns-x-y) + s~ V + e y ) - t(e~ x + e""), (A.4) 

such that 

^_9/ A c? _«9/a r _ 9/a t _ dfx 

5 "aT' T -~W (A ' 5) 

An alternative proof [10] is based on the fact that £{u) and M(u, A) are symplectic leaves of 
the same Poisson algebra (j2.23j) . 

Relation (|A.2|) defines thus a canonical transformation from XxSs to YyTt. 

Let us take n copies Wg^ of Ws decorating the variables XxYySsTt with the indices 

j = 1, ... , n and impose on them n matrix relations obtained from (1A.2|) by adding subscript j to 

(i) 

all variables. We obtain then a Lagrangian manifold Q^ n = <S>™ =1 £?4 in the 8n-dimensional sym- 
plectic space Wsn = ©?=iWg with the symplectic form = and the corresponding 

canonical transformation with the generating function ^ fx(yj,tj',Xj,Sj). 

i=i 

Let us also introduce 4 additional variables To, to an d S n +\, serving as coordinates in the 
4-dimensional symplectic space W4 with the symplectic form to>4 = dS n+ \ A ds n+ i — dTo A dto. 
The relations 



2( ai + (hXtp) _ 2(As n+ i - a r 



J ~~ 1 I a + 2 ' ~~ a TT2 \ A -^J 



define then a 2-dimensional Lagrangian submanifold Q2 C W4 characterised by the generat- 
ing function ip = ip x °\to) + <^ n+1 ^(s n+ i) with tp^ and V 9 !™" 1 "^ defined by (j4.5bj) and (|4.5c|) . 
respectively: 

T _ dip x dfx (A 7) 

dto os n+ i 

We end up with the (8n + 4)-dimensional symplectic space Wsn+4 = WW + W4, symplectic 
form Wgn+4 = ^8n+^4 ; and the (4n+2)-dimensional Lagrangian submanifold Qi n+ 2 = ^4n X^C 
Wsn+4 defined by the generating function 

n 

Fx = <p x °\t ) + <p[ n+1) (s n+l ) + £ fx{yj,t j]Xj , Sj ). (A.8) 

j'=i 

The final step is to impose 2n + 2 constraints 

T j = S j+u tj = s j+ i, j = 0, . . . , n, (A.9) 

which define a subspace W6n+2 C Wsn+4 °f dimension (8n + 4) — (2n + 2) = 6n + 2 and respective 
2n-dimensional submanifold C/2n = Qau+2 H Wq u+ 2- 

Constraints (|A.9|) allow to eliminate the variables Tt. The space W6n+2 splits then into the 
direct sum Wq u+ 2 = Vau + W211+2 of the space W4 n with coordinates XjXjYjyj (j = 1, . . . , re) 



14 



V. Kuznetsov and E. Sklyanin 



and W2n+2 with coordinates SjSj (j = 1, . . . ,n + 1). Using (|A.9|) we obtain that dTj A dtj — 
dSj + \ A dsj + i = and therefore the symplectic form (J8n+4 restricted on W§ n+ 2 

n 

U8n+4\ W6n+2 = ( dX J A dx J ~ dY 3 A Vj) > ( A - 10 ) 

i=i 

degenerates: it vanishes on W2 n +2 and remains nondegenerate on Vn n . In fact, on V^ n the form 
coincides with the standard symplectic form f|4.1j) . 

w 8n+4|y 4n = ^4n- (A. 11) 

After the elimination of the variables Tt from equations (|A.3[) and (|A.6|) . the resulting set 
of equations defining the submanifold Qm = Qin+2 H Wen+2 C Wen+2 coincides with equa- 
tions (IXT3]) and ([3T4]) defining the BT. 

As we have seen in Sectional the variables SjSj can also be eliminated leaving a 2n dimen- 
sional submanifold I^n C V& n coinciding with the graph of the BT discussed in Section 14.21 
By construction, I^n is the projection of Qm from Wqu+2 onto Vi n parallel to W2n+2- Further- 
more, T2n is Lagrangian since W8n+4 vanishes on £/4 n +2> therefore on Q2 n = Gin+2 H WQ n+ 2, and 
therefore on I^n- The canonicity of the BT is thus established geometrically, without tedious 
calculations. 

The same argument as in |10| shows that the generating function $>\ of the Lagrangian 
submanifold r 2n is obtained by setting tj = Sj + i in f|A.8|) . which produces formula (|4.4p . 



B Proof of spectrality 

Here we provide the proof of formulae ([4. 8ft for the eigenvalues of L(X). For the proof we use 
an observation from [10] and show that the eigenvectors of L(X) are given by null-vectors of 
Mi(±X,X). 

After setting u = — X in (|3.12p the matrix Mj becomes a projector 

x)= (-2X + s jSj ^-2A.A = /-2A + Si 5A(l s s ) ^ 



Sj SjSj J \ Sj 

with the null-vector 

Mj(-X,X)aj = 0. (B.2) 

Let us set u = — X in the matrix equality (13. ip and apply it to the vector a\. By ()B.2p . the 
right-hand side gives 0. Therefore, L{—X)a\ should be proportional to the same null-vector a\ 
of Mf(— A, A), and <ti is an eigenvector of L(— A). 

To find the corresponding eigenvalue A, use the factorised expression (|2.ip of L(— A) and 
apply it to o"i. Using (|2.3p we obtain 

^(-A; Y^yjVj = -Sje~ Vi aj + i, (B.3) 

hence 

n 

T(u;y,y)ffi = a n+1 n(-^e-w). (B.4) 

i=i 

From (]3.5p and A3. 12H we obtain 

A > = (7a; -Vsf u - ~\ + SjSj) • < B ' 5 > 
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hence 



the corresponding null- vector being 

A direct calculation using (|2.3|) and (|3.13d|) yields 

^(A: )). . i •' •rr / (B.8) 

and, consequently, 

n 

T t (X;Y,y)a n+1 = o x \\[s^). (B.9) 

i=i 

From (|2.4p we get, respectively, the identities 



I a 2 + PiX 2 

K + (-X)a n+ i = -((3 n + sl +1 )a n+1 , K-(-\)ai = 2 ^ + ^2 g i- 



(B.10) 



Using the above formulae we are able to move <ti through all the factors constituting L{— A) 
and obtain the equality 

L(-A;y,i/)<7i=Affi, (B.ll) 

where A is given by (|4.8ap . Note that A is an eigenvalue of L(A) as well since A (A) = A (—A). 
The second eigenvalue A (|4.8bj) of L(X) is obtained from 

AA = det L(A) = d(X) = {a 2 n + l3 n X 2 ){a\ + f5iX 2 ), (B.12) 

see fl233]). 

C Derivation of the dual Lax matrix 

To construct the 'big' Lax operator C(v) from the 'small' one L{u) we use the technique developed 
for the periodic the periodic Toda lattice [10l[19], with the necessary corrections to accommodate 
the boundary conditions. 

Let 61 be an eigenvector of L{u) with the eigenvalue v: 

L(u)e 1 = v6 1 , *i=fcY ( ai ) 

Reading off the factors constituting the product L(u), see (12. ip . (|2,2|) . define recursively the 
vectors 6j 



'3 



j = l,...2n + 2, (C.2) 



by the relations 

0j +1 = £(u;Xj,Xj)0j, j = l,...,n, (C.3a) 
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e n+2 = K + ( u )e n+1 , (c.3b) 

n+ j+3=e t (-u;X n -j,x n -.j)6 n +j+2, 3 = 0, ...,n-l, (C.3c) 
and close the circuit with the equation 

v9 1 = K_{u)e 2n+2 , (C.3d) 

which is equivalent to (jC.ip . 

A recursive elimination of ipj results in the equations 



(C.4a) 



\ V V J V V 

utfj = <p j+l - Xjipj + e x i- x i- Vn-i, j = 2,...,n (C.4b) 

u<p n+ i = ip n+ 2 + a n e~ Xn (p n , (C.4c) 
mp n+2 = -¥n+3 + X n ip n+2 + (a n e~ Xn + (3 n e~ 2xn )(p n+ i 

- (3 n e- 2x "X nVn + A.e^'-'-Vi, (C4d) 

uipj = -ip j+ i + X 2n+ 2-j<Pj ~ e^- 3_ ^- 4 ^-i, 3 = n + 3, . . . , 2n + 1, (C.4e) 

w<^2n+2 = u^i - aie x V2n+i- (C.4f) 

In order to simplify the 6-terms relations (|C.4aP and (|C,4dp and make the matrix C(v) more 
symmetric we perform an additional reversible change of variables 

<pi = <?i + 7r-e 2xi ^2n+i, (C.5a) 
2v 

( Pj = <Pji J = 2, ...,n + l, (C.5b) 

0n 

f n+2 = <p n +2 + -Ye~ 2Xn (f n , (C.5c) 

ifj = —<Pji 3 = n + 3, • • • , 2n + 2. (C.5d) 
The resulting equations for ipj read 

uft = £ 2 - Xift - f — e* 1 + ^e 2 * 1 ) £ 2n+2 - ^ +x ^2n, (C.6a) 
\ f 2v J 2v 

u& = <P3- X 2 y 2 + z X2 - Xl $i + ^e 2xi &n+i, (C.6b) 

wpj = tpj+i - Xj&j + e x ^~ x ^- 1 cp n - 1 , j = 3,... ,n (C.6c) 

n^ n+ i = + (a n e~ Xn + ^-eT 2x A £ n , (C.6d) 



Wp n+2 = -Vn+3 + X n Vn+2 ~ Ue Xji + ye - ye X " ^Vn-l, (C.6e) 

u^j = -<pj+i + X 2n+2 -j^ - e^- 3 "^- 4 ^_i, j = n + 4, . . . , 2n + 1, (C.6g) 

w^2n+2 = - ^aie Xl + ye 2:Cl ^ c^n+l- (C.6h) 

Introducing the vector O with 2n+2 components tpj, j = 1, . . . , 2n+2 we can rewrite relations 
(|C.6P in the matrix form 

C(v)S = u£(v)0 (C.7) 

with the matrix C(v) given by ipnj) . It follows from (107]) that u is an eigenvalue of 
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The rest of the formulae of Section [5] are obtained by a straitforward calculation not much 
different from the periodic case . 
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